® 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(tj) Numerode publication: 



0 348 247 

A1 



® 



DEMANDE DE BREVET EUROPEEN 



® Numerode depot: 89400854.9 
@ Date de depot: 24.03.89 



® Int. CI. 4 : G 01 B 11/00 

G 01 B 11/24, G 05 B 19/42, 
G 05 B 19/403, B 23 Q 35/128 



Priorite: 25.03.88 FR 8803991 

Date de publication de la demande: 
27.12.89 Bulletin 89/52 



Etats contractants designes: 

AT BE CH DE ES FR GB GR IT LI LU NL SE 



© Demandeur: KREON INGENIERIE MARKETING 
Route des Usines 
F-65300 Lannemezan (FR) 

@ Inventeur: Everaere, Michel 
Lotissement du Bois Pinas 
F-65300 Lannemezan (FR) 

@ Mandataire: Rinuy, Santarelli 
14, avenue de la Grande Armee 
F-75017 Paris (FR) 



Le titre de {'invention a ete modifie (Directives relatives a I'examen pratique a I'OEB, A-ill, 7.3) 

@ Procede de determination et de reconstitutlon des coordonnees spatiales de chacun des points d un ensemble de points 
echantillonnant une surface tridimensionnelle et reconstruction de cette surface a partir desdites coordonnees. 



(g) Le procede est notamment applicable a la video sculpture. 
On enregistre au moyen de cameras (13a, 13b) et de 
magnetoscopes (14a, 14b) une s6rie d'images de la surface 
tridimensionnelle eclairee par un faisceau plan de lumiere laser 
g6nere par un laser (12). Le procede comporte une phase de 
passage des coordonnees des points images enregistres par le 
systeme de vision aux coordonnees spatiales mise en oeuvre 
par ie module (17), une phase de reconstitution des points 
eventuetlement rnal enregistres par le systeme de vision, mise 
en oeuvre par le module (16) et une phase de determination de 
I'axe d'un outil (18) sculptant I'image tridimensionnelle dans 
une matiere telle que du bois ou de Taluminium, cette derniere 
phase etant mise en oeuvre par le module (19). 
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Description 

Procede de determination et de reconstitution des coordonnees spatiales de chacun des points d'un 
ensemble de points echantillonnant une surface tridimensionnelle, et procede de realisation d'une image 
tridimensionnelle de cette surface a partir desdites coordonnees. 

5 La presente invention concerne, d'une maniere generate, le probleme de la realisation d'une image 
trimensionnelle d'une surface quelconque a partir d'images bi-dimensionnelles de ladite surface. Elle trouve 
une application importante, mais non exclusive, dans le domaine de la "video sculpture". 

L'invention vise notamment un procede de determination des coordonnees spatiales de chacun d'un 
ensemble de points appartenant a ladite surface trimensionnelle et echantillonnant cette derniere, a partir 
10 desdites images preenregistrees. Grace aux donnees procurees ainsi, l'invention permet une commande 
optimale d'une machine-outil a commande numerique pour usiner ladite image tridimensionnelle dans un 
rnateriau, dans le cadre de son application a la "video sculpture". 

On connait un tel procede permettant, dans i'art anterieur, de realiser I'image tridimensionnelle d'un buste 
a partir d'une serie d'images video. Le sujet est assis sur un fauteuil tournant et son buste est eclaire au moyen 
15 d'un laser fournissant un faisceau plat iamellaire. L'intersection du faisceau plat avec le buste du sujet donne 
une courbe (ou profil) que Ton echantillonne, les coordonnees spatiales des points source des echantillons 
etant determinees au moyen d'un procede de type connu, tel que celui enseigne dans le brevet francais 
81.24418. Le fauteuil sur lequel est assis ie sujet effectue une rotation pas a pas et a chaque fois une courbe 
est ainsi numerisee. En une rotation complete, on obtient ainsi un nombre predetermine de courbes (par 
20 exemple 512) echantillonnees chacune en un nombre predetermine de points (par exemple 512). Le buste est 
ainsi nurnerise sous la forme des coordonnees X-Y-Z de chacun des points memorisees (au nombre de 
262.144 dans i'exemple ci-dessus). 

L'etape suivante consiste a traiter les points numerises ainsi obtenus de facon a generer des ordres 
permettant de commander la machine-outil, de telle sorte que cette derniere scuipte dans une rnatiere solide 
25 telle que de Taluminium ou du bois i'image tridimensionnelle du buste. 

La presente invention apporte notamment un perfectionnement au procede decrit ci-dessus. 
Le procede connu jusqu'a present presente, en effet, un certain nombre d'inconvenients, notamment dans 
la phase de traitement des informations et de generation des signaux electriques susceptibles de commander 
la machine-outil, cette phase etant communement appelee dans les milieux specialises, phase de CFAO 
30 (conception et fabrication assistees par ordinateur). 

Ainsi, le procede d'eclairement du buste n'est pas parfait en ce sens que la peau ou les cheveux du sujet 
reflechissent la iumiere emise par le laser de facon variable, ce qui peut se traduire par des points sureclaires 
ou des points sous-eclaires. Dans certaines circonstances, la peau ou les cheveux du sujet peuvent presenter 
une faible albebo, ce qui signifie que le point eclaire ne reflechit pas de Iumiere. Dans un tel cas, la trace 
35 enregistree par la camera video comporte des absences de points creant des ruptures de continuity de la 
trace. On considere alors qu'il existe des points "manquants" qui auraient du combler les vides ou les ruptures 
de con tinuite. Ces points "manquants" correspondent aux points de buste presentant un faible albedo. 

A Tinverse, des points de fort albedo peuvent reflechir la Iumiere de facon tres forte et diffractent la Iumiere 
sur d'autres parties du buste, la trace enregistree n'etant pas continue et comportant en fait divers points 
40 exterieurs a la trace normale. II s'agit dans un tel cas de points que Ton qualifle "d'aberrants". 

Dans Tun et I'autre cas, I'usinage, sous le controle de la CFAO, manque de precision, ce qui genere des 
erreurs prejudiciables pour la qualite du produit fini. 

Une autre serie d'inconvenients est egalement due au fait que le procede utilise pour determiner les 
coordonnees spatiales a partir des coordonnees des points image enregistres par la video, fait appel a des 
45 formules trigonometriques telles que celles decrites dans le brevet 81.24418 mentionne plus haut. II s'ensuit 
un manque de precision. II est done souhaitable de disposer d'un procede de determination des coordonnees 
spatiales des points considered, plus precis et plus fiable. 

Par ailleurs, la Demanderesse a constate qu'il etait souhaitable d'orienter I'outil de telle sorte qu'il attaque la 
rnatiere sous un angle correspondant a la normale a la surface devant etre realisee, ce qui permet d'obtenir un 
50 meilleur fini et un meilleur rendu de la piece. L'art anterieur ne prevoit pas une telle hypothese puisque la 
machine qui est utilisee n'est commandee que selon trois degres de liberie, ce qui ne permet que de 
positionner de facon correcte la fraise attaquant la rnatiere, mais ne permet pas de positionner Tangle 
d'attaque de cette fraise par rapport a la rnatiere. II est done souhaitable de disposer d'un procede permettant 
de commander et de controler autant que possible Tangle d'attaque de i'outil par rapport a la rnatiere. 
55 Enfin, la Demanderesse a constate que dans certaines zones telles que les yeux ou la bouche, si Tangle 
d'attaque de la fraise est normal a la surface a usiner, il risque de se produire une collision entre Tarbre 
porte-outil et la rnatiere usinee. Par consequent, dans ces hypotheses, il faut pouvoir rectifier, de facon 
automatique dans toute la mesure du possible, la position de Tarbre porte-outil de telle sorte que I'outil 
attaque la rnatiere sous un angle favorable mais que Ton evite, autant que possible, les collisions entre Tarbre 
60 porte-outil et la rnatiere. A cet egard, il est done souhaitable de disposer dans la CFAO d'une phase 
permettant de determiner la direction optimale de Tarbre porte-outil. 
La presente invention resout Tensemble de ces problemes. 

Elle vise notamment en effect un procede de determination de coordonnees spatiales de chacun des points 
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d'un ensemble de points source appartenant a une surface comportant une etape d'acquisition de donnees 
dans laquelle on 6ciaire avec un faisceau plan de lumiere la surface, on enregistre avec au moins une camera 
au moins une image de la trace laissee par le faisceau plan sur la surface, la trace laissee par le faisceau plan 
de lumiere sur la surface etant dite "profil source", rimage de cette trace dans un plan image de la camera 
etant dite "profil image", on depiace relativement le faisceau plan et la surface et on enregistre avec ladite 
camera a nouveau au moins une image de la trace laissee par le faisceau plan sur la surface, on recommence 
les operations de deplacement relatif, d'eclairage et d'enregistrement de sorte que I'on enregistre ainsi une 
pluralite de profils image constituant un echantillonnage de rimage de I'ensemble de ia surface, le precede 
comportant une etape de traitement au cours de laquelle, notamment, on considere les coordonnees 
lignes-colonne (Col.Lig) d'un ensemble de points image echantillonnant chaque image profil dans le plan 
image de ladite camera, procede caracterise en ce que, mettant en oeuvre les considerations suivant 
lesquelles : 

- dans le plan du faisceau lumineux, on choisit un axe directeur et on affecte des coordonnees dites 
"spatiales" (R,H), a chaque point dudit plan, R etant la distance entre le point considere et I'axe directeur, H 
etant I'ordonnee de ce point sur cet axe ; 

- on prend en compte les formules suivantes de transformation des coordonnees (Col,Lig) de chacun des 
points du plan image de ia camera en coordonnees spatiales (R, H) du point source correspondant : 

a 1 Lig+a 2 Col+a 3 a 4 Lig+a 5 Col+a 6 

(1) R = = ^ H « 

a 7 Lig+a 8 Col+i a 7 Lig+a 8 Col+l 

les coefficients a t etant dits "coefficients de transformation" et etant determines pour 1'ensemble des points 
du plan image de ladite camera, a partir de coefficients dits "d'interpolation" (ftj) eux-memes calcules lors 
d'une phase prealable de calibration du plan image de la camera au cours de laquelle : 
. on enregistre rimage de series juxtaposees de quatre points source de reference situes dans le plan de la 
lumiere et dont les coordonnees spatiales (R, H) dans ce plan sont connues ; on releve les coordonnees 
images (Col.Lig) de rimage de chacun de ces points pour chaque serie de quatre points source de reference 
et on resout la formule de transformation (1 ), ce qui donne huit valeurs de coefficients de transformation oti par 
serie de quatre points de reference ; on considere le barycentre de chacune des series de quatre points de 
reference et on affecte a ce dernier les coefficients de transformation ai ainsi determines ; on considere des 
zones d'interpolation Tk de forme triangulaire delimitees par des segments de droite joignant lesdits 
barycentres de series voisines, ces barycentres constituant les sommets desdites zones triangulaires T k ; on 
utilise la formule d'interpolation suivante pour chaque coefficient de transformation ar. 

(2) <* - Pii Col + p i2 Lig 4- p e 

. on determine vingt-quatre coefficients d'interpolation (pij) par zone triangulaire d'interpolation (Tk) en 
resolvant la formule d'interpolation (2) pour chacun des coefficients de transformation (ai) affecte a un 
sommet de zone triangulaire d'interpolation Tk ; 

ce procede consiste en outre, dans ladite phase de traitement, pour determiner les coordonnees spatiales 
(R,H) d'un point source dont on connait les coordonnees images (Col,Lig) de son image, a effectuer les 
operations suivantes : 

- rechercher, dans une operation de recherche, la zone d'interpolation (Tk) a laquelle le point appartient, ainsi 
que les coefficients d'interpolation associes a cette zone. 

- calculer, dans une operation de calcul, les coordonnees spatiales (R, H) par application concomitante 
desdites formules de transformation (1) et d'interpolation (2). 

Selon un autre aspect de I'invention, il est mis en oeuvre une phase d'homogeneisation au cours de laquelle, 
pour chacun des points traites (P,, n ), on considere un nombre predetermine de points voisins, on calcule les 
coordonnees moyennes de ces points, si le point traite est absent, ou si Tune au moins de ses coordonnees 
differe de la coordonnee moyenne correspondante d'un ecart superieur a une valeur predeterminee (E), on 
affecte au point traite lesdites coordonnees moyennes. 

Ainsi, grace a ces dispositions, on dispose d'une methode fiable et precise pour determiner et/ou 
reconstituer les coordonnees spatiales d'un ensemble de points a partir de leurs coordonnees dans le plan 
image de la camera. 

De plus, grace a la phase d'homogeneisation, on pallie les inconvenients lies aux points de faible ou de trop 
forte albedo sur la surface a reproduire. 

La presente invention vise egaiement un procede de video sculpture, notamment du type tel que 
succinctement rappele ci-dessus, pouvant, avantageusement, mettre en oeuvre le procede de determination 
et/ou de reconstitution des coordonnees spatiales enonce ci-dessus, notamment caracterise en ce qu'H met 
en oeuvre une machine-outil adaptee a commander un outil selon cinq degres de liberie, a savoir trois degres 
de liberie perrnettant le positionnement de cet outil dans l'espace, compte tenu desdites coordonnees 
spatiales et deux degr6s de liberte determinant un angle d'attaque de I'outil de telle sorte que I'outil soit 
generalement dirige selon la normale a la surface devant etre realisee et en ce que Ton procede a une 
detection des parties concaves de la surface, pour corriger Tangle d'attaque de I'outil pour eviter autant que 
possible une interaction entre I'outil et la matiere usinee, au cours de laquelle : 
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- on considere la trace laissee par ladite surface dans un plan perpendiculaire a un axe directeur et passant par 
le point en cours de traitement (Pi, n ), cette trace etant dite secante (SO, 

- on considere egalement un profii source 0 n auquel appartient le point traite (Ri, H„ 0 n ), 

- on detecte les parties concaves du profll en effectuant les operations suivantes : 

5 . on considere un vecteur normal ( DN* en, m, dn^ enj+i) a chacun des deux points voisins (Ri-i,HM,0n) et 
(R i+1 ,Hi + i,e n ), du point traite (Ri.Ht.On) 

. on calcule le produit vectoriel desdits vecteurs normaux ('DlT en,M, D*f enj+i) au profii aux points 
voisins du point traite ; on determine ainsi un coefficient de convexite QD, associe au point en traitement (Ri, 

Hi, 0n), 

10 . sHe coefficient de convexite est negatif, le profii source est concave au point traite, et on calcule la somme 
(__e5 en, 0 desdits vecteurs normaux au profii, aux deux points voisins du point traite, ( DN*" en,M, 

DN 6n,i+l), ^ 

. on memorise ce vecteur somme ( ED en.i), 

. on considere aiors I'ensemble des points traites appartenant a un mime profii (6 n ) ; on associe a chacun de 
15 ces points soit le vecteur normal au profii au point considere, lorsque le coefficient de convexite QD 
considere est positif, soit le vecteur somme, lorsque le coefficient de convexite (QD) est negatif, 
. lorsque I'ensemble des points appartenant a un meme profii est traite, on recommence un_ nombre 
predetermine de fois le traitement du profii en utilisant les vecteurs effectivement associes ( Dlt e n , i, 
ElJ en.O au point traite (R if H, 6n). 
20 - on procede a un traitement analogue pour chacune des traces secantes (Si) et on associe a chacun des 
points traites (Rj, Hi, On) un vecteur correspondant soit a la normale a ladite trace secante Si au point traite 
SN enj, soit a un vecteur somme, ES en.O des vecteurs effectivement associes aux points voisins (Ri, Hi, 
0 n -i) (Ri, Hi, 8 n +i) au point traite (Ri, Hi, On). 

et, pour chaque point, on dirige I'outil selon ia direction des vecteurs effectivement associes audit point, dans 
25 la place de son profii (0 n ) et dans celui de la trace secante (S-,) a laquelle il appartient. 

Grace a ces dispositions, on contrdle de facon optimale I'outil, tout en evitant, de facon automatique le plus 
souvent, toute interference entre I'outil et la matrere. 

Les caracteristiques et avantages ressortiront d'ailleurs de la description qui va suivre en reference aux 
dessins annexes surlesquels : 
30 - la figure 1 est une vue schematique des moyens mettant en oeuvre un procede de video sculpture 

conforme a ! 'invention, 

- la figure 2 est une vue schematique illustrant la position relative du laser et des cameras illustres en 
figure 1, 

- les figures 2a et 2b iliustrent une variante de realisation, 

35 - la figure 3 est une vue schematique d'une surface tridimensionnelie a reproduire, 

- la figure 4 illustre, de fapon schematique, un profii image tel qu'enregistre par chacune des cameras 
de la figure 1, 

- la figure 5 illustre le profii source correspondant au profii image de la figure 4, 

- la figure 6 est une vue schematique d'une mire de calibration du systeme de vision, 

40 - ia figure 7 est une vue schematique d'une autre mire de calibration permettant une calibration plus 

precise de chacune des cameras, 

- la figure 8 est un schema fonctionnel des principaux moyens mis en oeuvre dans le mode de 
realisation decrit et represente, 

- la figure 9 est un organigramme simplifie de transformation realisant I'une des etapes illustrees en 
45 figure 8, 

- La figure 10 est une illustration de cinq profils source du sujet illustre en figure 3, 

- les figures 11 a 13 iliustrent un profii source presentant une partie concave et les etapes permettant 
de corriger la direction de Toutil pour eviter une interaction outil-matiere et, 

- les figures 14 et 15 iliustrent une trace "secante" presentant une partie concave. 

50 On va tout d'abord brievement decrire un mode de realisation du procede selon I'invention et les moyens le 
mettant en oeuvre pour avoir une vue generale. On decrira ensuite plus en detail certaines phases du procede. 

1. Generalites 

la figure 1 illustre de facon generale des moyens mettant en oeuvre un procede de video sculpture 
55 comportant le procede de realisation d'une image tridimensionnelie selon la presente invention. 

Le procede de video sculpture illustre comporte, de fapon classique, trois etapes principales mises en 
oeuvre par des ensembles de moyens : 

- I'ensemble A permet ['acquisition de ('information, 

- I'ensemble B permet le traitement de I'information, 

60 - I'ensemble C est une utilisation de I'information, en I'espece usinage de I'objet copie au moyen d'une 

machine-outil 18. 

L'acquisition de I'information se fait en placant I'objet a copier sur une table 10 ou le sujet dont on desire 
modeler le buste sur un siege 10'. La table 10 et le siege 10' sont rotatifs et entraines a cet effet par un moteur 
pas a pas non represente. Le sujet ou I'objet sont eclaires au moyen d'un laser 12. Dans le mode de realisation 
65 choisi et represente, un seul laser eclaire le sujet ou I'objet. Deux cameras 13a, 13b filment I'objet ou le sujet 
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alors que la table 10 ou le siege 10' effectue un tour complet. Les pas angulaires de la table ou du siege sont ici 
de 1° en sorte que deux rnagnetoscopes 14a, 14b respective ment raccordes aux cameras 13a, 13b 
enregistrent ainsi chacun 360 images p de I'objet ou du sujet eclaire par le laser 12 (voir figure 3). 

Sur la figure 3, on a illustre un sujet 20 et deux traces p iaissees par le laser 12. Le double systeme 
d'acquisition des donnees 13a, 14a et 13b, 14b enregistre ainsi une pluralite de traces ou "profils" p realisant, 5 
de ce fait, un echantillonnage de I'image de I'ensembie de la surface du buste 20 du sujet. 

Avant de poursuivre la presente description, il est donne les definitions suivantes : 

- "point source" : point appartenant a la surface exterieure de I'objet ou du sujet que Ton veut reproduire, 

- "profii source" : trace p laissee par ie laser sur le buste 20 ou I'objet a reproduire, 

- "point image" : image d'un point source enregistre par la camera 13a et/ou la camera 13b, dans un "plan 10 
image" appartenant a la camera consideree. Le point image est caracterise par ses coordonnees (Col,Lig) 
dans le plan image ainsi qu'il ressortira du reste de la description, 

- "profii image" : image du profii source dans le plan image de la camera consideree. 

La figure 2 illustre schematiquement la disposition relative des deux cameras 13a, 13b et du laser 12. 
L'ensemble est vu du dessus, la tete du sujet portant la reference 20, 1'axe 120 du laser 12 et celui (130a, 130b) 15 
de chacune des cameras 13a, 13b se coupe sur un axe vertical Z correspondant a l'axe de rotation du siege 10' 
sur lequel est assis le sujet. L'axe Z est dit, dans la suite de la description, "axe directeur". On observe que 
I'axe Z debouche de ia tete 20 sensiblement au milieu de cette derniere. Les axes 130a, 130b, sont inclines 
d'un angle 8 de 30° par rapport a l'axe 120 du laser 12. Ce dernier emet un faisceau de lumiere plan qui laisse 
une trace, egalement appelee profii et referencee sur la figure 2 en p, sur ia tete et sur le buste du sujet. 20 

On choisit Tune des deux cameras, ia camera 13a en I'espece, comme camera directrice. En regie generale, 
cette derniere voit en permanence et dans sa totalite la trace du faisceau laser sur la tete et le buste du sujet. 
Cependant, dans certaines zones, la trace p est partiellement cachee, et ne peut etre vue de la camera 
directrice 13a (voir par exemple la figure 2 : au niveau du nez, la trace p est, pour une portion de son etendue 
tout au moins, masquee). Aussi, dans un tel cas, on utilisera les donnees issues de la camera secondaire 13b, 25 
qui, s'agissant de I'exemple illustre en figure 2 tout au moins, voit la portion masquee de la trace p. Plus loin, 
dans la description, il sera explique comment un fichier directeur Fla, associe a la camera directrice 13a, est 
complete par les donnees issues de la camera secondaire 13b. 

II convient ici d'observer que le mode d'acquisition des donnees est particulierement bien adapte a une 
surface presentant une forme variable sur les 360° de sa peripherie. 30 

Cependant, une telle disposition n'est pas imperative. En effet, au lieu de placer I'objet a reproduire sur un 
plateau tournant, on pourra, par exemple, le disposer sur un plateau anime d'un mouvement de translation de 
telle sorte que les traces seront alors paralleles entre elles et paralleles a un axe directeur solidaire de I'objet, 
cet axe se deplacant parallelement a lui-meme. Ce pourra &tre le cas lorsque I'objet a reproduire est par 
exemple une carte en relief dont la representation mathematique est de la forme : z = f(x, y) et dont seuie une 35 
face est a reproduire. En revanche, lorsque I'objet doit etre reproduit dans son ensemble, on utilisera le 
procede du plateau tournant decrit plus haut. 

A propos du mode de realisation illustre, on notera par ailieurs que, avantageusement, compte tenu de la 
forme de I'objet, ou de la morphologie du sujet, il est possible de parametrer les deplacements angulaires pour 
obtenir plus d'images sur les zones de I'objet du sujet dont la surface est sujette a de brusques variations (par 40 
exemple le nez du sujet) et moins d'images dans les zones ou la surface est plus reguliere (par exemple, la 
nuque du sujet). 

Dans I'exemple illustre, les donnees sont enregistrees sur un premier ensemble de rnagnetoscopes 14a, 
14b, un autre ensemble de rnagnetoscopes 14'a, 14'b etant utilise pour le traitement de I'information. Ceci 
permet de decentraliser la phase de traitement par rapport a celle d'acquisition. 45 

En figure 2a, il est illustre une variante de realisation de I'ensembie d'acquisition de donnees mettant en 
oeuvre une seule camera 13c et deux lasers 12a, 12b. Cette figure illustre schematiquement la disposition 
relative des deux lasers et de la camera. L'ensemble est vu du dessus, la tete du sujet portant aussi la 
reference 20. L'axe 130c de la camera et celui (120a, 120b) de chacun des lasers 12a, 12b se coupent sur un 
axe vertical Z correspondant a i'axe de rotation du siege 10' sur lequel est assis le sujet. Les axes 120a, 120b 50 
sont inclines de 30° par rapport a I'axe 130c de la camera. Les lasers 12a, 12b emettent chacun un faisceau de 
lumiere. Ce faisceau laisse sur la tete du sujet et son buste une trace, referencee sur la figure 2a en p (trace 
laissee par le laser de gauche 12a) et en p' (trace laissee par ie laser de droite 12b). 

En regie generale, la camera voit en permanence les deux traces des lasers sur la tete de la personne. Ces 
deux traces sont completement separees car elles ne peuvent avoir de point commun que sur l'axe de rotation 55 
Z. Les moyens electroniques ne prennent en consideration que I'une des traces, en I'espece la trace p', 
laissee par le laser de droite 12b. Cependant, comme prec6demment, dans certaines zones, la trace de droite 
sera cachee (zone masquee par ie nez du sujet, par exemple). Aussi, dans un tel cas, la camera 13c ne peut 
voir cette trace ou tout du moins la portion de celle-ci qui est masquee par le nez. Dans ce cas, les moyens 
electroniques prennent en compte la trace de gauche vue par la camera lorsque la tete a pivote de 60° comme 60 
illustre en figure 2b. En effet, la trace p laissee par le laser de gauche 12a (de la figure 2b) est identique a celle 
(p') laissee par le laser de droite 12b de la figure 2a. Grace a cette disposition, on peut ainsi compenser et 
retrouver les zones masquees par le nez ou les oreilles par exemple. 

Cette variante de realisation presente pour avantages d'alleger sensiblement le dispositif d'acquisition des 
donnees puisque Ton fait I'economie d'une deuxierne camera et de deux rnagnetoscopes. Cependant, cette 65 
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variante de realisation presente, par rapport au mode prefere de realisation choisi, deux inconv6nients : d'une 
part, la mise en oeuvre de la prise en compte des 60° de decalage complique senslblement le traitement et, 
d'autre part, les contraintes mecaniques sont telles que le profil p' (figure 2a) ne coincide pas parfaitement 
avec le profil p (figure 2b) une fois que la tete a pivote de 60° : en effet, le sujet est incapable de rester 
5 parfaitement immobile pendant ce temps. C'est pourquoi, en definitive, le Demanderesse a prefere mettre en 
oeuvre le mode de realisation decrit a I'appui des figures 1 et 2. Cependant, la variante illustree a I'appui des 
figures 2a et 2b demontre qu'il existe divers procedes permettant d'enregistrer I'ensemble de la surface a 
reproduire, y compris les parties cachees. Aussi, la presente invention englobe-t-elle a cet egard toutes 
variantes susceptibles d'arriver a ce resultat. 

10 Le signal analogique issu des cameras 13a, 13b enregistre dans les magnetoscopes 14a, 14b et delivre par 
les magnetoscopes 14'a, 14'b se presente sous la forme d'une serie de courbes (les profils p). Les courbes 
sont sensiblement continues et presentent des points tres eclaires et des points faiblement eclaires. 

Chacune de ces courbes est numerisee dans les moyens electroniques 15, de facon connue en soi, au 
moyen d'un convertisseur analogique-numerique qui delivre les coordonnees image (ligne, colonne), d'une 

15 suite d'echantillons representatifs de chacune des courbes. Le brevet francais 81.24418 decrit un mode 
prefere de realisation permettant de realiser cette numerisation, et par consequent le convertisseur 15 ne sera 
pas decrit plus en detail ici. 

Cependant, la Demanderesse a constate que pour des raisons d'ordre technique et economique, Tobjet ou 
le sujet a reproduire etant dispose sur un support rotatif, les coordonnees les plus facilement exploitables et 

20 les plus economiques du point de vue taille memoire a utiliser sont les coordonnees dites "cylindriques" (R ,H, 
9), R correspondant a la distance entre un point de la surface et Taxe de Z tandis que H correspond a 
I'ordonnee le long de I'axe des Z et 9 correspond a la coordonnee angulaire du profil source (intersection d'un 
axe X' appartenant a un referentiel solidaire du buste 20, figure 3, avec I'axe X dans le plan XY fixe). 
La coordonnee 9 correspond en fait au pas (ou a un nombre de pas predetermine et parametrable) du 

25 moteur commandant la rotation du siege 10'. Une impulsion est enregistree dans la bande son des 
magnetoscopes 14a, 14b a chaque pas du moteur ; comme par ailleurs, la bande son est synchronisee avec la 
prise d'image, on connait la position du moteur (et par consequent Tangle 9) a chaque prise d'image. 

Les moyens de traitement B de I'invention comportent, outre les magnetoscopes 14'a, 14'b et le 
convertisseur 15 des moyens 17 procedant a une reconstitution des coordonnees cylindriques des points 

30 source a partir des coordonnees image de chacun des echantillons composant les profils. Le precede de 
reconstitution mis en oeuvre par les moyens de transformation 17 est decrit plus loin. 

En sortie des moyens de transformation 17, on dispose dinformations portant sur les coordonnees 
cylindriques d'un nombre predetermine de points representatifs de la surface tridimensionnelfe, ici celle du 
buste 20 a reproduire. 

35 Les moyens de traitement de I'information B comportent un ensemble de CFAO compose de deux modules 

16 et 19 raccorde aux moyens de transformation 17. Cet ensemble traite les coordonnees issues des moyens 

17 et permet la commande d'une machine-outil 18. Cet ensemble met en oeuvre, notamment par des moyens 
informatiques, un procede decrit en detail plus loin. 

Les moyens de transformation 17, et I'ensemble CFAO 16, 19 sont, dans le mode de realisation prefere 
40 presentement decrit, constitues par un micro-ordinateur 21 realisant les fonctions des moyens 17, 16 et 19 
decrites plus loin et contrdlant, de surcroTt, le convertisseur 15. Le micro-ordinateur 21 est ici un ordinateur 
IBM, modele commercialise sous la reference AT3. 

Enfin, la machine-outil 18 qui, ici, est du type NC-TRIAX-R6-2-TWIST, fabriquee par la Societe italienne CMS, 
permet de sculpter dans un materiau tel que du bois, de la resine ou de I 'aluminium, une forme reproduisant 
45 sensiblement celle du buste du sujet assis dans le siege 10' ou celle de I'objet a reproduire dispose sur la table 
10. 

En resume, dans le mode de realisation presentement decrit, le procede traitant des donnees issues du 
convertisseur 15 comporte trois phases : 

- une phase de passage des coordonnees image des points image enregistres aux coordonnees cylindriques 
50 de leurs points source. C'est la fonction des moyens de transformation 17, 

- une phase d'homogeneisation des points, (fonction du module 16), 

- une phase de determination de la normale a la surface et de correction de cette normale, pour determiner la 
direction de I'outil (fonction du module 19). 

Dans la description qui va suivre, on va decrire en detail chacune de ces trois phases. On commencera par 
55 expliquer tout d'abord la structure des donnees fournies par le systeme de vision constitue ici par les 
cameras, les faisceaux laser, les magnetoscopes et le convertisseur. 

2. Structure des donnees en sortie du convertisseur 15 
Chacune des 360 images enregistrees correspond a un profil source, caracterise par i'angle 9 n forme entre 
60 I'axe X' et I'axe des X dans le plan X-Y. 

On notera dans la suite de la description que 6 designe tout a la fois la coordonnee angulaire d'un profil 
source, le plan contenant ce profil source caracterise par sa coordonnee angulaire 9, et le profil image 
correspondant au profil source 9. La variable 9 est indicee notamment selon n, n correspondant au numero 
d'ordre du profil source ou image considere. 
65 Chaque camera video 13a, 13b est une camera par exemple de type THOMSON TAV-1033-NEWICON 
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delivrant 316 lignes par image et 1024 pixels par ligne. Aussi, les profiis image 6 n sont-ils en fait constitues par 
une succession de points caracterises par leur coordonnee dans i'image video enregistree, ces dernieres 
correspondant en fait a celles du pixel. En effet, pour un pixel donne, on connaTt ia ligne sur laquelle il est 
dispose (coordonnee Lig variant de 0 a 315) et la position du pixel dans cette ligne (coordonnee Col variant 
de 0 a 1023). " 5 

Un point quelconque du profil image 8 n est done defini par ses coordonnees : (Col, Lig, 8 n ). 

En figure 4, on a illustre de facon schematique, un profil image 6 n tel qu'enregistre par la camera video. On a 
schematise par une croix ie pixel eclaire a I'occasion de la prise d'images. Chaque image est constitute par 
316 lignes ordonnees de 0 a 315, 1024 pixels etant susceptibles d'etre eclaires sur une ligne. On a indique, a 
titre indicatif, pour le profil 0 n illustre, le numero du pixel eclaire au-dessus de la croix schematisant ce pixel. 10 
Ainsi, dans I'exemple illustre en ligne 1, le pixel n° 480 est eclaire, en ligne 3, le pixel n° 500 est eclaire, etc. 

Pour chaque point du profil 0 n , on dispose done de ses coordonnees image (Col, Lig, 9 n ). A titre d'exemple, 
celles du point situe sur la sixieme ligne de la figure 4 sont : (502,6, 8 n ). 

Les moyens de transformation 15 mettent en oeuvre un fichier image Fl dans lequel les coordonnees image 
des points constituants les divers profiis sont stockees de ia maniere suivante : 15 

- le fichier image contient sequentiellement les ordonnees des pixels du profil, sur une ligne consideree, 

- chaque coordonnee est codee sur un mot de 16 bits. 
Le code d'un point image se fait de la facon suivante : 

TABLEAU I 20 



N° du bit 


Commentaire 


Bit 15 


0 separateur de profil 




1 point du profil en cours 


Bit 14 


0 pas de point dans la ligne en cours 




1 presence d'un point dans la ligne 




en cours 


Bit 13-Bit 10 


0 ces bits ne servent pas 


Bits 9 - 0 


Position du pixel eclaire dans la 




ligne codee sur 10 bits 



25 



30 



A chaque profil image correspond un ensemble sequentiel de 316 mots. Les profiis sont separes par le mot: 
0 0 0000 0000000000. Le numero du profil en cours de lecture est inscrit dans un compteur. 35 
Le tableau II ci-dessous illustre a propos du profil 6 n la structure du fichier Fl de la figure 4. 

TABLEAU II 

0 0 0000 0000000000 Separateur de 40 

profil 

Profil 0n 1 0 0000 0000000000 Lig 0, pas de 

point 

1 1 0000 01111000000 Lig 1, col 480 

1 0 0000 0000000000 Lig 2, pas de 45 
point 

1 1 0000 0111110100 Lig 3, col 500 
1 1 0000 0111101010 Lig 4, col 490 
1 1 0000 1011010000 Lig 5, col 720 

1 1 0000 0111110110 Lig 6, coi 502 50 
1 1 0000 0111010110 Lig 7, col 470 
1 1 0000 0111100110 Lig 308, col 486 
1 1 0000 0111110101 Lig 309, col 501 

1 0 0000 0000000000 Lig 310, pas de 55 
point 

1 1 0000 0111100101 Lig311,col485 
1 1 0000 0011111010 Lig 312, col 250 
1 1 0000 0111100110 Lig313,col486 

1 1 0000 0111110100 Lig 314, col 500 60 
1 1 0000 0111101111 Lig315,col495 
0 0 0000 0000000000 Separateur de 
profil 

Profil Gn + i 1 1 0000 0111100011 Lig 0, col 493 ^ 
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En se referant a la figure 4 et au tableau II, on observe que chaque ligne peut comporter soit : 

- 1 point dans une zone mediane a, 

- 1 point dans des zones laterales b, 

- aucun point dans ces zones. II s'agit en fait de points insuffisamment eclaires pour ies raisons expliquees 
5 ci-apres. 

Ainsi, en lignes 2 et 310, aucun pixel n'est eclaire. En revanche, en lignes 5 et 312, des pixels respectivement 
720 et 250 se trouvent situes dans Ies zones laterales b. Ces points sont dits "aberrants". 

En effet, le systeme de vision utiiisi est tel qu'il est possible que certains points de Tobjet ou du sujet ne 
refiechissent pas suffisamment la lumiere laser et, meme, au contraire I'absorbent. Ces points sont des points 
10 dits "de faible albedo". Dans un te! cas, la lumiere reffechie etant insuffisante, le pixel concerne ne se trouve 
pas eclaire. C'est pourquoi en lignes 2 et 310, dans I'exempie illustre en figure 4, aucun pixel ne se trouve 
eclaire. 

En revanche, il est possible que certains pixels se trouvent eclaires, ces pixels ne correspondant pas a des 
points du sujet appartenant au prof il eclaire : ce sont ici Ies points "aberrants" des lignes 5 et 312 (par exemple 
15 pour cause de reflet de la lumiere laser dans Ies cheveux du sujet). 

On decrit plus loin une caracteristique de invention permettant de reconstituer Ies points manquants et de 
detecter Ies points aberrants. 

En definitive, on observera qu'en sortie du convertisseur 15 et en entree des moyens de transformation 17, 
on dispose d'un fichier recapitufant de facon sequentielle Ies coordonnees (Col, Lig, 9) d'ensembles de 316 
20 points appartenant a 360 profils. 

3. Reconstitution des coordonnees spatiales du point source de chacun des points image enregistres 

Le probieme pose est celui de permettre le passage des coordonnees image d'un point (Col, Lig, 9) aux 
coordonnees spatiales de son point source, ces coordonnees spatiales etant en t'espece choisies de type 
25 cylindrique (R, H, 9). 

Ce probieme n'est pas simple car il s'agit de passer d'un systeme de coordonnees dans un plan 
particulierement reduit, celui de la cellule de la camera, a des coordonnees spatiales, en trois dimensions. 

Une premiere methode consiste a determiner, de facon mathematique, pour chaque point du plan image Ies 
coordonnees cylindriques (R, H, 0) du point source correspondant, etant entendu que Tangle 0 est celui defini 
30 plus haut, ces coordonnees etant mernorisees dans un fichier et la transformation coordonnees image (Col, 
Lig, 0) en coordonnees cylindriques (R, H et 9) se faisant par simple lecture du fichier. 

Cette methode, bien qu'elle entre dans le cadre general de la presente invention, n'est pas satisfaisante 
techniquement et economiquement au moment du depot de la presente demande, car elle exige une taille 
memoire particulierement importante. II est cependant possible qu'elle soit techniquement et economique- 
35 ment envisageable dans un proche avenir et que la presente invention englobe une telle methode. 

La Demanderesse a par consequent utilise et ameliore une technique permettant de transformer Ies 
coordonnees image en coordonnees cylindriques utilisant une taille memoire relativement restreinte. 

En fait, s'agissant des coordonnees cylindriques (R, H, 9), le probieme se pose essentieilement pour la 
determination des coordonnees R et H, Tangle 0 etant Tangle entre Taxe des X' et Taxe des X (figure 2), cet 
40 angle etant donne par la position du moteur pas a pas animant le siege 10'. 

Le modele mathematique suivant est utilise : 

^Lig+o^Col+o^ a 4 Lig+a 5 Col+a 6 

45 " R - H = ■ 

a 7 Lig+a 8 Col+l a 7 Lig+a 8 Col+l 



Le probieme qui se pose done id est done de determiner Ies coefficients de transformation oti a appliquer 
50 aux coordonnees ligne-colonne de tous Ies points de Tespace susceptibles d'etre filmes par Ies cameras 13a, 
13b, pour obtenir Ies coordonnees cylindriques correspondantes pour chacun de ces points. 

A cet effet, on procede a une etape de " calibration " de chacune des cameras 13a, b, in situ. C'est Tobjet de 
la phase de calibration decrite ci-apres au paragraphe 3.1. 
Une fois ces coefficients de transformation determines comme i! sera decrit ci-apres, on peut proceder 
55 grace aux equations indiquees sous (1) ci-dessus a la transformation proprement dite des coordonnees 
ligne-colonne des points delivres par le convertisseur 15 pour chacun des systemes de vision 13a, 13b. On 
genere ainsi deux fichiers dits "spatiaux" FSa et FSb (Cf. paragraphe 3.2. ci-apres). 

3.1. Phase de calibration 

60 La methode de base utilisee pour la calibration est la suivante : 

On utilise une mire de calibration telle qu'illustree en figure 6, comportant une plaque rectangulaire 41 de 
hauteur H et de largeur R dans Texemple considere R - H = 50 cm) et d'epaisseur e (1 cm ici). La plaque 41 
est munie de quatre piges constitutes par des cylindres metalliques (referencees A, B, C, D dans Texemple) 
de rayon r et de hauteur h predetermines (dans Texemple considere r = 1 mm et h = 2 cm). On connait la 

65 position spatiale des piges A - D et, notamment, leurs coordonnees (RA, HA), (RB, HB), (RC, HC) et (RD, HD). 
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Sur la figure 6, on a illustre I'axe des Z, le long duquel les coordonnees H sont mesurees et un axe R, 
perpendiculaire a I'axe des Z, le long duquel les coordonnees R peuvent etre mesurees. On a egalement 
illustre, sur cette figure, les coordonnees (RA, HA) de la pige A. 

On filme la mire ainsi constitute au moyen de chacune des cameras 13. On a illustre sur cette figure, sous la 
reference 42, I'axe de vision de la cam6ra 13a et sous la reference 43, le plan lamellaire de lumiere issu du laser. 5 
Ce plan lamellaire est parallele a la plaque 41 , tandis que I'axe de vision 42 de la camera 13 s'incline de 30° par 
rapport au plan lamellaire de lumiere et intercepte la plaque 41 au barycentre 0 des quatre piges A - D. 

Le systeme de vision enregistre done les coordonnees image correspondant aux quatre piges : (ColA, 
LigA), (ColB, LigB), (ColC, LigC) et (ColD, LigD). 

A partir des equations (1) ci-dessus, il est possible de former un systeme de huit equations a huit inconnues 10 
(les a.) et comme parmi les quatre points A - D, il n'existe aucun sous-ensemble de trois points alignes, le 
systeme d'equations admet une solution a la fois unique et non aberrante. 

Le modele mathematique (1) ci-dessus utilise avec les ai ainsi determines donne a I'evidence une solution 
exacte pour les quatre points A-D. Cependant, a cause des non-linearites optiques du systeme de vision, la 
precision des coordonnees RH pour les autres points se degrade. On a observe que I'erreur est en fait 15 
minimale au barycentre des autres points A-D. 

Pour minimiser cette erreur, selon une caracteristique de cet aspect de I'invention, on procede de le facon 
suivante. On utilise une mire de calibration se presentant sous la forme d'une plaque carree 53 dont les 
dimensions correspondant a la plus grande dimension de I'objet ou du sujet que Ton souhaite reproduire en 
trois dimensions (voir figure 7). Dans I'exemple considere, la mire presente une hauteur et une largeur de 1 m. 20 
On dispose une pluralite de piges 51 en ligne et en coionne de facon a definir une pluralite de carres 52. Dans 
le mode de realisation choisi et represents on divise ainsi la plaque 53 de calibration en 25 carres sur chacun 
des coins desquels une pige de calibration 51 est disposee. Chaque pige est caracterisee par ses 
coordonnees ligne - coionne (i, j). Ainsi, la pige (i, j) est la pige de la ieme ligne et de la jeme coionne. 

On met en place la mire de calibration 50 comme montre en figure 6. On a illustre sur la figure 7 un axe des Z 25 
et un axe R respectivement paraileles a deux cotes perpendiculaires de la mire de calibration 51. En 
connaissant les coordonnees (R,H) des piges 51, on peut determiner, pour chacun des barycentres 54 de 
chacun des carres 52, les coefficients a\ afferents au carre considere, en appliquant les equations (1) 
ci-dessus. On affecte alors a chacun des barycentres 54 les coefficients cti correspondants. 

On considere ensuite les triangles T k dont les sommets sont constitues par les barycentres 54 de carres 30 
voisins. Ces triangles T k sont illustres en figure 7 en traits mixtes. Les triangles T k sont ici au nombre de 32 et 
sont numerates (k variant de 1 a 32). 

Pour assurer une continuity du modele mathematique utilise, il est considere que les coefficients en sont des 
polynomes de degre 1 en Col, Lig sur chacun des triangles T k . Ainsi, les coefficients a, sont de la forme : 
(2) ai - PmCoI + pisLig + Pi3 35 

Pour chacun des triangles T k , on connait les coefficients ai permettant de determiner les coordonnees R, H 

de chacun de ses sommets. 

On dispose done, pour chaque triangle T k> et pour chaque a\ afferent a un sommet, d'un systeme de trois 
equations a trois inconnues (les ft) que Ton peut resoudre pour determiner ces inconnues. 

En resume, la methode est done la suivante : 40 

- pour chacun des barycentres 54, on determine les coefficients ct< permettant de determiner les coordonnees 
spatiales (R, H), compte tenu des coordonnees image (Col, Lig) enregistrees par le systeme de vision, et 
compte tenu du fait que i'on connait les coordonnees spatiales (R, H) des piges 51 visionnees par le systeme 
de vision, 

- a partir des coefficients a\ affectes a chacun des barycentres, on determine les coefficients ft pour les 45 
triangles T k considered 

Par consequent, on dispose ainsi d'une methode d 1 interpolation permettant de faire vaner continuellement 
les coefficients on, compte tenu de la position des points disposes a I'interieur d'un triangle T k donne et, 
notamment de leur rapprochement ou de leur eloignement d'un sommet considere. 

Avec chacune des cameras 13a, 13b, on filme la mire 50 une fois pour toutes, et on procede aux calculs des 50 
coefficients pi, afferents a chacune des cameras. 

Une fois que les coefficients pij ont ete determines pour les deux cameras, celles-ci sont alors "calibrees" et 
on peut transformer les coordonnees (Col, Lig) issues du convertisseur 15 en coordonnees spatiales (R,H). 

3.2. Constitution des fichiers associ6s a chaque camera & 

Comme explique plus haut, chaque camera 13a, 13b filme une serie de profils source, chacun des profils 
image resultant etant echantillonne par le convertisseur 15 de telle sorte qu'en sortie de ce convertisseur, on 
dispose des coordonnees (Col.Lig) d'un ensemble d'echantillons permettant de reconstituer chacun des 
profiis consid6res. 

La figure 8 est un schema fonctionnel des principaux moyens mis ici en oeuvre. La fonction des cameras 60 
13a, 13b et ceile de leurs magnetoscopes associes 14a, 14'a, 14b, 14'b y ont ete represents de facon 
schematique. Sur cette figure, on a egalement represent6 schematiquement les fonctions principaies 
effectuees par le convertisseur 15 et les moyens transformation 17. 

Le convertisseur 15 comporte essentiellement des moyens de conversion analogique-numerique 
fonctionnant selon le principe decrit dans le brevet francais 81.24418 cite plus haut. Ces moyens de 65 



9 



EP 0 348 247 At 



conversion portent la reference 151. Les coordonnees (Col.Lig) pour chacun des points des profils image 
delivrees par chacune des cameras sont en fait enregistrees dans deux fichiers image Fla, Fib. Le fichier Fla 
est associe au systerne de vision comportant la camera 13a, le rnagnetoscope 14a et ie magnetoscope 14'a, 
tandis que le fichier Fib est associe au deuxieme systerne de vision (camera 13b, magnetoscopes 14b et 14'b). 
5 Ainsi, seules les coordonnees des points en provenance de ce second systerne de vision sont enregistrees 
dans le fichier Fib, tandis que seules les coordonnees des points issus du premier systerne de vision sont 
enregistrees dans le fichier Fla. 

Les moyens de transformation 17 comportent un fichier de transformation 171a (171b) associe a chacune 
des cameras 13a, 13b. Chaque fichier 171a (171b) contient un nombre d'enregistrements egai au nombre de 
10 triangles T k considere (ici 32). Chaque enregistrement comporte : 

- le numero du triangle Tk considere, 

- les coordonnees de chacun des sommets de ce triangle (ban/centres 54), 

- les coefficients p y associes a ce triangle pour la camera consideree 13a(13b) (obtenus lors de la calibration 
de ces cameras). Comme expose plus haut, la calibration de chaque camera est faite in situ prealablement a 

15 toute exploitation commerciale du precede. Sur la figure 8, la calibration des cameras, permettant de 
constituer les fichiers 171a, 171b, est schematisee en C a et C b ). 

Les moyens de transformation 17 comportent egaiement, d'une part un fichier spatial FSa constitue des 
coordonnees cyiindriques obtenues apres transformation des coordonnees image contenues dans ie fichier 
image Fla et, d'autre part, un fichier FSb constitue des coordonnees cyiindriques obtenues apres 
20 transformation des coordonnees image dans ie fichier image Fib. 

A partir des coordonnees image d'un point quelconque, contenues dans Tun des fichiers Fla ou Fib, on 
caicule, au moyen d une methode iterative, les coordonnees spatiales R, H correspondant a ce point. 

En figure 8, I'etape iterative de transformation a ete schematisee en 172a, pour les coordonnees issues du 
fichier Fla, et en 172b pour les coordonnees issues du fichier Fib. La figure 9 represente un organigramme 
25 simplifie de transformation, realisant I'etape 172a ou 172b. On ne decrira que I'etape 172a, I'etape 172b etant 
identique. 

On traite chacun des points contenus dans !e fichier Fla de facon sequentielle. Le traitement d'un point 
determine, ayant des coordonnees (Lig, Col) comporte une lecture dans le fichier Fla, et une inscription en 
memoire, comme schematise en 180. L' etape de transformation decrite vise a calculer les coordonnees (R, H) 

30 correspondant aux coordonnees (Lig, Col). Pour ce faire, lors d'une premiere etape d'initialisation 181, on 
initialise au depart de chaque traitement d'un point image determine, en mettant la memoire qui doit contenir 
les coordonnees du point source (R,H) a zero. Le processus iteratif est alors le suivant : en 182, lors d'une 
seconde etape d'initialisation, on initialise deux variables intermediaires (R Ql H D ) en mettant ces variables 
intermediaires egales au contenu de la memoire (R,H). La variable K, qui correspond au numero du triangle 

35 est, elleaussi, initialisee a zero. En 183, on incremente cette variable k. En 184, on vient lire dans le fichier 171a 
les coordonnees du triangle Tk, ainsi que le jeu du coefficient (3y associe au triangle T k considere. 

En 185, on procede a un test d'appartenance pour determiner si le point ayant pour coordonnees les 
variables intermediaires (R c , H c ) appartient ou non au triangle T k . Pour ce faire, on fait des tests de 
comparaisons successives de coordonnees selon une methode de recherche classique bien connue de 

40 rhomme de Tart qui sait, connaissant les coordonnees d'un point dans un plan determine, a quelle region de 
ce plan appartient ce point. Si le test 185 est negatif, on revient en 183 ou Ton incremente la variable k. Lorsque 
le test est positif, e'est-a-dire que i'on a determine le triangle T k auquel appartient le point ayant les 
coordonnees intermediaires (R D , Ho), on passe a une etape de transformation 186, 
En 186, on affecte aux coordonnees (R, H) contenues en memoire de nouvelles valeurs correspondant a la 

45 transfor'mee des coordonnees (Col,Lig), au moyen des equations (1) et (2) ecrites plus haut, permettant de 
calculer a partir des coefficients lus en 184 les coordonnees (R, H) considerees. 

En 187, on procede au test consistant a determiner si les coordonnees (R, H) sont respectivement egales 
ou non aux coordonnees intermediaires (R 0 , H G ). Dans la negative, on revient a I'etape 182 ; si le test 187 est 
positif, cela signifie que les coordonnees (R, H) en memoire sont effectivement celles du point source du point 

50 image (Lig, Col) traite. 

On inscrit ces coordonnees (R, H) dans le fichier spatial FSa. 

En resume, au cours de I'etape iterative de transformation, il est effectue les operations suivantes : 

- une operation de recherche au cours de laquelle on recherche la zone d'interpolation T k a laquelfe le point 
source appartient, ainsi que les coefficients d'interpolation Pi] associes a cette zone. On observera qu'au lieu 

55 d'effectuer cette recherche dans I'espace source, comme dans le mode de realisation decrit, il est possible, 
sans sortir du cadre de I'invention, de I'effectuer dans le plan image de la camera 13a ou 13b, 

- une operation de calcul au cours de laquelle on caicule les coordonnees spatiales (R,H) par application 
concomitante des formules de transformation (1) et d'interpolation (2). 

On constitue de meme un fichier spatial FSb a partir des coordonnees contenues dans le fichier Fib. 
60 il convient d'observer que pour un meme point spatial, les coordonnees cyiindriques determinees a partir 
des coordonnees image dans le fichier Fla et dans le fichier Fib devraient etre les memes. En pratique, elles 
seront identiques aux erreurs optiques et d'arrondi pres. 

On observera que le procede de transformation decrit ci-dessus utilise peu de memoire puisqu'en I'espece, 
il suffit de memoriser par triangle T k : 
65 - son numero (une information), 
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- les coordonn6es spatiales de ses sommets (six informations), 

- les pij assoctes (24 informations), 

soit trente et une informations par triangle. II y a done un total de 992 informations a memoriser, ce qui est 
reiativement faible a comparer aux informations qui devraient etre memorises s'il fallait memoriser dans une 
seule memoire les coordonnees transformees pour chacun des 1024 pixels de chacune des 316 lignes des 5 
cameras 13, soit au total 323.584 informations). 

Les coordonnees spatiales contenues dans les fichiers FSa et FSb sont fusionnees en 173 en un seul fichier 
FS. Cette fusion permet de completer avec les coordonnees des points de I'un des fichiers les coordonnees 
absentes dans I'autre. 

C'est a ce niveau du procede qu'est solutionne le probleme des points qui n'ont pu etre fiimes par les deux 10 
cameras (parties cach6es) : les parties cachees pour une camera ont pu etre filmees par I'autre (voir figure 2). 
En procedant a la fusion des fichiers FSa et FSb, on realise un fichier unique FS dans lequel en theorie 
('ensemble de la surface du buste est 6chantillonne (sous reserve des points manquants ou aberrants dans les 
deux fichiers pour les raisons mentionnees plus haut). 

La fusion s'opSre en utilisant I'un des fichiers, le fichier FSa, comme fichier pilote et en completant les 15 
coordonn6es manquantes dans ce fichier par celles existant Sventuellement dans le fichier FSb. 

4. Homogeneisation des points d'un profil 

Comme explique plus haut, a I'appui de la figure 4, chaque profil B n peut : 

- soit presenter des absences de points, 20 

- soit comporter des points aberrants. 

Le but de I'etape presentement decrite, qui est dite "d'homogeneisation", est d'eliminer les points aberrants 
et de reconstituer les coordonnees des points manquants. Cette etape comporte une operation dite 
n d'int6gration", une operation de "filtrage" et une operation dite de "reconstitution". 

25 

4.1 Operation dlntegration 

En figure 10, il est illustre cinq profils source (8 n -2, G„-i, 9 n , 8„+i, 9 n +2). On observe que ces profils 
concernent cinq positions successives du siege 10'. lis ont ete illustres dans le referentiel associe au buste 20 
(X', Y', Z). 

Le profil source en cours de traitement est le profil 0 n . Pour chacun des points traites : 30 
Pi, n - (Ri, Hi,9n), on effectue les operations suivantes : 

- on determine tout d'abord les points voisins du point (Ri, Hi, 9 n ) consider^. A cet egard, on considere les 
2p + 1 profils aliant de 9 n - P a e n+P . Dans I'exemple illustre en figure 10, p = 2 et par consequent, on considere 
les cinq profils allant de 9 n -2 a 9 n +2. 

Sur chaque profil source 0j ainsi considere, on selectionne les 2k + 1 points allant de (Rh<, Hj-k, 0j) a (Ri+k, 35 
Hi +k , 9]). Dans I'exemple considere, en figure 10, on a choisi pour valeur de k : k = 2. 

Cela determine ainsi (2p + 1) x (2k + 1) points voisins de (Ri, Hi, 8 n ) y compris ce point en cours de 
traitement. 

Dans le mode de realisation represente, il y a done un total theorique de 25 points voisins au point Pi, n en 
cours de traitement, y compris ce dernier. _ 40 

On calcule ensuite la moyenne R des ordonnees en R et la moyenne H des ordonnees en Z de chacun de 
ces 25 points. 

Dans i'hypothese ou des points sont manquants, la moyenne est effectuee en diminuant le diviseur commun 
(25 dans un cas normal), d'autant d'unites que de points sont manquants. 
Si le point "traite" Pi, n n'est pas manquant (bit 14 a 1 dans le tableau I), on effectue I'operation de fiitrage. 45 

4.2 Operation de filtrage 

Le filtrage a pour but de permettre Telimination des points aberrants. Pour chacun des points traites, on 
calcule la valeur absolue I Ri - R I . 

Si cette valeur est superieure a une valeur d'erreur donnee E, parametrable, on elimine le point (Ri, Hi, 8 n ). 50 

La valeur donnee E est determinee par Texperience, compte tenu de I'objet que Ton cherche a reproduire. 
S'agissant par exemple du visage d'une personne, pour certaines zones, on sait que le risque d'obtenir des 
points aberrants est particulierement eleve et que ces points aberrants risquent d'etre disposes a une 
distance reiativement importante du profil. Dans un tel cas, on choisira pour valeur E, une valeur de quelques 
centimetres, 3 cm par exemple. 55 

En revanche, dans d'autres cas, les erreurs et les points aberrants seront peu nombreux et en tous cas 
reiativement proches du profil. Dans un tel cas, on choisira comme valeur E, une valeur faible de I'ordre de 
quelques millimetres, 3 mm par exemple. 

On observera en plus, selon une caracteristique de I'invention, que ladite valeur donnee E peut etre 
parametrable compte tenu du profil analyse. Ainsi, dans certaines zones du visage, telles que le nez, on 60 
choisira une valeur E reiativement importante, puisque le profil est susceptible d'avoir d'importantes variations 
de sa coordonnee R, tandis que, pour les profils concernant la nuque, on choisira une valeur E reiativement 
faible, puisque ie profil n'a que des variations de sa coordonnee R de faible amplitude. 

65 
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4 3. Operation de reconstitution v tii . 

Dans le cas oO le point (R, H,. Bh) "trait6" est manquant (bit 14 a 0, Cf. tableau I), ou s H est filtre comme 
expose ci-dessus, on effectue I'operationdite de reconstitution au cours de laquelle on affecte au po.nt (R,, H,, 

9 n ), les valeurs moyennes calculees R, H, lors de I'etape d'initialisation. 

5 La figure 5 est une vue correspondant a la figure 4, dans laquelle les points representes sont les points 
source en coordonnees cylindriques, et dans laquelle les points aberrants et manquants ont ete reconstitues, 

comme expose ci-dessus (points P r ). ,j„„„a D o h= 

Le fichier FS est complete au fur et a mesure du traitement des points (Ri, Hi 6„) par les coordonnees de 

tous les points P r reconstitues. ... ■ 

10 On notera que I'etape d'homogeneisation d'un profil 9 n qui vient d'etre decrrte met en oeuvre les 
coordonnees cylindriques, apres reconstitution des coordonnees cylindriques du point source a partir des 

coordonnees de son image. . _ 

Dans un autre mode de realisation de la presente invention, II est tout a fart envisageable d mtervertir les 
etapes de reconstitution et d'homogeneisation et de proceder a I'etape d'homogeneisation en coordonnees 

15 image puis de proceder ensuite a une transformation ou k la reconstitution des coordonnees cylindriques de 
I'ensemble des points source a partir des points image d'origine ou des points "image" reconstitues. On 
observera en effet que, d'une maniere generate, la phase d'homogeneisation est une phase au cours de 
laquelle pour chaque point traite (que ce dernier sort considere dans I'espace source ou dans le plan image), 
on considere un nombre predetermine de points voisins, on calcule les coordonnees moyennes de ces points, 

20 si le point traite est absent, ou si I'une au moins de ses coordonnees differe de la coordonnee moyenne 
correspondante d'un ecart superieur a une valeur predetermine* (E), on affecte au point traite lesdites 
coordonnees moyennes. 

5. Dete rmination de la direction de I'arbre porte-outil _ 
25 L'outil utilise se presente e ssentiellement sous la forme d'une fraise disposee en bout d un arbre 
appartenant a un robot susceptible de commander la position de cet arbre dans I'espace. 

On concoit que la position de la fraise doit etre determine* au moins partrois coordonnees. En I espece, les 
coordonnees (R, H, 9) seront utilisees. 
Mais il faut de surcroTt, que I'arbre portant la fraise soit dispose selon une direction sensiblement normale 
30 a la surface a reproduire, de telle sorte que la matiere soit attaquee sous un angle favorable 

Pour ce faire, le robot doit etre capable de determiner non seulement la position de la fraise selon les trois 
coordonnees cylindriques, mais en plus la direction generale de I'arbre (deux degree de liberte 

SU [Ce porte-outil sera oriente selon une direction dite normale a la surface si deux conditions se trouvent 

35 X -Tux\t part, II faut que la projection de I'arbre dans le plan 9„ auquel appartient le point traite: Pi, n : (Ri, Hi, 9„) 
soit normale au profil auquel appartient ce point, 

- et d'autre part, que la projection de I'arbre dans un plan perpendiculaire secant Si, passant par le point traite 
Pi „ soit egalement normale a la trace de la surface a reproduire dans ce plan secant. Dans la suite, on appellera 
40 "trace se-cante" la trace laissee dans ie plan Sipar la surface a reproduire. 

En figure 1 1 , on a illustre sous la reference DN^i la projection dans le plan 9„ de la normale a la surface au 
point P, n , tandis qu'en figure 15, on a represents sous la reference SN Q „,i la project.on de la normale a cette 
surface au point Pi >n dans le plan Su . , „ . „ ... . „ aU ^* 

Ainsi, en plus des trois degres de liberty necessaires pour commander la posit.on de axe de I outil le > robot 
45 qui commande cette position, ainsi que la position de I'arbre, doit etre capa ble de ^I tZ ^ 

Hberte supplemental. Les deux coordonnees supplemental que le robot dort ma.triser sont des 
rotations. L'une autour de I'axe des Z (perpendiculaire au plan S,). I'autre autour d'un axe perpendiculaire a 

i'axe des Z, I'axe Y en I'espece. , ,-. „„ o»iic 

Le robot et l'outil utilises sont ici du type NC-TRIAX-R6-2-TWIST, fabrique par la Societe ttalienne CMS 
Un probleme se pose lorsque la surface a usiner comporte plusieurs parties concaves, I une de celles-c. 
etant representee en figures 12 et 14. En effet, dans une telle partie concave, il risque de se produ.re une 
collision outil-matiere. Ainsi, figure 1 2, au point critique Pc, la direction de la projectton de la normale i est D Or 
cette direction intercepte le profil non seulement au point Pc mais egalement au point Pin de telle sorte que si 
I'arbre de l'outil prenait une telle direction, il y aurait collision avec la matiere. 
55 Le meme probleme se pose d'ailleurs dans le plan Si, figure 14. „ moHan »H'i«iter 
Pour pallier cet inconv6nient, il faut determiner une direction dite de ''normale rectifiee permettant d eviter 
cette collision et, dans ce but, s'ecartant de la normale. Cette direction a ete illustree pour le point Pc en figure 
13 sous la reference "ifienj et en figure 14 (dans le plan SO sous la reference ES Q n.i. En consequence, 
I'etape de determination de la direction de l'outil se decompose en trois sous-etapes : 
60 - estimation de la normale N 9n .i, . . 

- detection des parties concaves dans le plan 9 r et calcul de la direction de la projection des normales 

rectifiees dans le plan 9 n , , , .. . 

- detection des parties concaves dans le plan Si et calcul de la direction de la projection des normales 

rectifiees dans le plan Si. 



SO 
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5.1. Estimation de la normale 

La determination de la normale du point traite ft.n va se faire, compte tenu des quatre points les plus 
proches P2-P5, figure 10. Pour determiner cette normale : on procede au calcul des produits vectoriels 
suivants 



L 
> 



Ni - Pt?2 A P t p 4 

N" 2 
N 
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3 ~ 
N 4 " = P t P5 

On norme les vecteurs Hi*- N 4 T ce qui donne les vecteurs 15 

NMi - "N*T 4 (N^i = — r- ) 

llNill 

On calcule la somme des quatre vecteurs normes: 

25 

NMl + N«2 + + *™4 

ce qui permet d'obtenir un vecteur de direction moyenne aux quatre vecteurs normes precedents. La direction 
du vecteur resultant est consideree comme etant la "normale" a la surface au point en cours de traitement Pi, n . 30 
Le vecteur resultant est ensuite norme, ce qui donne le vecteur onj. 

Le traitement ci-dessus s'effectue pour chacun des points (Ri, Hi, 8 n ) du fichier FS. _^ 
On cr6e un fichier FN comportant, pour chacun des points Pj, n , les coordonnees du vecteur associe N en.i. 

5.2 Determination des parties concaves dans 8 n et calcul de la projection de la normale rectifiee dans ce plan 35 
On considere les vecteurs et les coordonnees suivants : 

- DN* 9 n j : projection normee de n anj dans le plan 8 n , 

- RNt et HNi : les composantes de ce vecteur dans le plan 6 n (figure 11). 

On considere les deux points voisjns de (Ri, Hi.JM respective me nt (Rm, Hm, e n ) et (Ri+i, Hi+i, 0 n ) sur le 
profil e n et la projection normee DNq nj-i. et dNq n,i+i associee a ces points. 40 
On calcule le produit vectoriel suivant : 



(^n / i+l A DN^ enfi ^ 1 ) 45 

ce qui permet d'obtenir un coefficient de convexite QDi associe au point en traitement (Ri, Hi, 0 n ) : 
QDi - (RNi+1 * HN i-i) -{RNm * HNt+i) (3.1) 

Si le coefficient QDi est positif, cela signifie que le meridien est convexe au point traite. S'il est negatif , cela 
signifie que Se meridien est concave. K 50 

Si le meridien est convexe, on inscrit dans le fichier FS les coordonnees du vecteur dn on.i associe au 
point traite (R i} Hi, 8 n ). 

Si le coefficient QDi est negatif, on calcule le vecteur : 

DNVi-1 + DN*6n.i+1 (3.2) 

On norme ce vecteur, ce qui donne le vecteur ed en.i > 55 

On inscrit dans le fichier FS les coordonnees du vecteur ED en.u 
On traite tous les points Pj. n du fichier FS. 

On dispose a ce moment d'un fichier FS comportant, d'une part les coordonnees spatiales de chacun des 
points Pin et, d'autre part associees a chacun des points consideres, les coordonnees : 

- soit de la projection DN* en.i de la normale N** enj dans le plan 6 n , 60 

- soit du vecteur en.i qui est une estimation de la direction de la projection de la direction de I'arbre 
porte-outil dans le plan 9 n . pour minimiser le risque d'interference entre porte-outil et matiere, mentionne plus 
haut. ED*en,i est appele "normale rectifiee dans le plan 6 n ". 

Lorsque I'ensemble des points est traite, on recommence le traitement ci-dessus, en utilisant, dans les 
formules (3.1) et (3.2) ci-dessus, soit les coordonnees du vecteur DN^nj ou Etfenj associees a chacun 65 
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des points et on modifie, dans le fichier FS, les valeurs des vecteurs ED^n,!. 

Ce processus est recommence un nombre predetermine de fois (trois dans cet exemple). En realite, la 
direction ErTgn.*, estimee au point considere apres trois iterations, est telle que le risque de collision 
outil-matiere mentionne plus haut, du a la concavite du meridien au point considere, est minimise. 

5 

5.3 Determination des parties concaves dans Si et calcul de la projection de la normale rectifiee dans ce plan 

On rappelle tout d'abord que le plan S\ est un plan perpendiculaire aux divers plans 8 n et interceptant I'axe 
des Z au niveau de la coordonnee Hi, figures 14 et 15. Ce plan est parallele au plan XY illustre en figure 3. 
On considere egalement la trace secante laissee par la surface du buste 20 dans le plan Si : Si designant 
10 done soit le plan, soit la trace secante. 

On considere les vecteurs et les coordonnees suivantes : 
SNen.i : projection normee de n@ n.i dans le plan Sj, 
- RSj et YSj les composantes de ce vecteur dans le plan Sj. 

On considere un point traite (Ri, H i} 8 n ). On considere les deux points voisins de ce point respectivement (Ri, 
15 Hi 9 n -i) et (Ri, Hi, 9 n -n) sur la trace secante Si et la projection normee sSqVu et SN« n +i i associee a ces 
points. 9 y 

On calcule le produit vectoriel suivant : 



sn en-l,i S N en+i,i 

ce qui permet d'obtenir un coefficient QSi associe au point en traitement Ri, Hi % : 
QSj - (RSi+i * YSi + i) - (RSj+1 * YSm) (4.1) 
25 Si le coefficient QS, est positif , cela signifie que la trace secante est convexe au point traite. S'il est negatif, 
cela signifie que la trace est concave. 

Si la trace est convexe, on inscrit dans le fichier FS les coordonnees du vecteur SN% n i associe au point 
traite (R n , Hi, G n ). 

Si l e coefficient QSi est negatif, on calcule le vecteur 

30 SN%n-l.i + SN**9n+1,l (4.2) _ 

On nomme ce vecteur, ce qui donne le vecteur es*W 

On inscrit dans le fichier FS les coordonnees du vecteur ES%n,i. 

On traite tous les points P iin du fichier FS. 

On dispose a ce moment d'un fichier FS comportant, d'une part les coordonnees spatiales de chacun des 
35 points, et les coordonnees de la projection dans le plan 8 n de la normale reelle ou rectifiee, (respectivement 
DNgnj et EDQn.i, et, d'autre part, les coordonnees : 

- soit de la projection^ sN%n.i de la normale de N% n ,i dans le plan Si, 

- soit du vecteur ES^en.i qui est une estimation de la direction de la projection de la direction de I'arbre 
porte-outil dans le plan Si minimisant le risque d'interference entre le porte-outil et la matiere mentionnee plus 

40 haut. E^enj est appele "normale rectifiee dans le plan Si". 

Comme precedemment, ce processus est recommence un nombre predetermine de fois (trois dans cet 
exemple). Ainsi, la direction ESgnj estimee au point considere, apres trois iterations, est telle que le risque 
de collision outil-matiere mentionne plus haut est minimise. 

45 5.4 Determination finale de la direction de I'arbre porte-outil 

Lorsque I'etape definie au paragraphe 5.3 precedent est terminee, le fichier FS comporte finalement pour 
chaque point Pj, n : 

- les trois coordonnees du point Pi, n (Ri, H f , 8 n ), > 

- les deux coordonnees dans le plan 9 n du vecteur DN 0 ni ou du vecteur rectifie EDQ n>i . Ces deux 
50 coordonnees definissent la direction de I'outil dans !e plan 6 n . v 

- les deux coordonnees dans le plan Si du vecteur SNQn.i ou du vecteur rectifie ESQ nl . ces deux 
coordonnees definissent la direction de I'outil dans le plan Sj. 

Ce fichier est alors transfere dans la machineroutil 18. 

La machine-outil 18 sera done piiotee a partir des trois coordonnees (Ri, Hi 9 n ) et des deux directions 
55 susmentionnees. II est a noter que la machine-outil CMS utilisee ici est tout a fait adaptee pour etre piiotee de 
facon automatique des lors que les coordonnees ci-dessus lui sont delivrees sous un format equivalent et 
approprie bien connu de Thomme de Tart, notamment celui familier des post-processeurs. 

Grace a ces caracteristiques, I'outil attaquera la matiere selon Tangle le plus favorable tout en evitant, ou, a 
tout le moins, en minimisant les risques de collisions porte-outil - matiere mentionnes plus haut. 
60 II faut toutefois observer que dans certains cas (rares) de forte concavite, les coefficients QDj et QSi 
calcules apres trois iterations ne seront pas positifs, et que dans de tels cas, une collision outil-matiere est 
susceptible de se produire. La Demanderesse a constate que ces cas sont rares, et, en pratique, s'agissant de 
video sculpture d'un visage, inexistants. Cependant, pour la reproduction de volumes tres tourmentes, il est 
possible de prevoir une alarme avec intervention manuelle permettant de corriger la direction de I'arbre 
65 porte-outil en cas de tres forte concavite (cas d'une cavite par exemple). II n'en reste pas moins que, grace a 
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I'invention, ces cas d 5 intervention manuelle restent fortement minimises. 

L'ensemble des etapes de caicul decrites depuis le paragraphe 3 ci-dessus, de meme que, d'une fapon 
generate, la gestion de l'ensemble du procede et notamment le controle et la commande des fonctions et des 
processus decrits a I'appui des figures 8 et 9, ainsi que celles de la machine outil, sont assures par le 
micro-ordinateur 21. Ce contrdle et cette gestion sont de nature classique, quoique adaptee au procede qui 5 
vient d'etre decrit. C'est pourquoi, ils ne seront pas decrits plus en detail ici. 

Bien entendu, la presente invention n'est nullement limitee aux modes de realisation decrits et represents. 
Bien au contraire, en plus des variantes deja indiquees, elle englobe toute autre variante a la portee de 
I'homme de I'art 



Revendications 

1 . Procede de determination des coordonnees spatiales de chacun des points d'un ensemble de points 
source appartenant a une surface comportant une etape d'acquisition de donnees dans laquelle on 15 
eclaire, avec un faisceau plan de lumiere, la surface, on enregistre avec au moins une camera (13a, 13b) 
au moins une image de la trace (P,p,0n) laissee par le faisceau plan sur la surface (20), la trace laissee par 
ie faisceau plan de lumiere sur la surface etant dite "profil source", I'image de cette trace dans un plan 
image de la camera etant dite "profil image", on deplace relativement le faisceau plan et la surface et on 
enregistre avec ladite camera a nouveau au moins une image de la trace laissee par le faisceau plan sur la 20 
surface, on recommence les operations de deplacement relatif, d'eclairage et d'enregistrement en sorte 
que Ton enregistre ainsi une pluralite de profils image constituant un echantillonnage de I'image de 
l'ensemble de la surface, le procede comportant une etape de traitement au cours de laquelle, 
notamment, on considere les coordonnees lignes-colonne (Col, Lig) d'un ensemble de points image 
echantillonnant chaque profil image dans le plan image de ladite camera, procede caracterise en ce que, 25 
mettant en oeuvre les considerations suivant lesquelles : 

- dans le plan du faisceau lumineux, on choisit un axe directeur (Z) et on affecte des coordonnees dites 
"spatiales" (R,H) a chaque point dudit plan, R etant la distance entre le point considere et I'axe directeur, 
H etant I'ordonnee de ce point sur cet axe : 

- on prend en compte les formules suivantes de transformation des coordonnees (Col, Lig) de chacun 30 
des points du plan image de la camera en coordonnees spatiales (R, H) du point source correspondant : 

Lig+a 2 Col+a 3 a 4 Iiig+a5Col+a 6 



a 7 Lig+a 8 Col+l otyliig+a^Col+l 



les coefficients ai etant dits "coefficients de transformation" et etant determines pour l'ensemble des 
points du plan image de ladite camera, a partir de coefficients dits "d'interpolation" (Py), eux-memes 
calcules lors d'une phase prealable de calibration du plan image de la camera au cours de laquelle : 40 
. on enregistre I'image de series juxtaposees de quatre points source de reference (51) situes dans le 
plan de lumiere et dont les coordonnees spatiales (R, H) dans ce plan sont connues ; on releve les 
coordonnees images (Col, Lig) de I'image de chacun de ces points pour chaque serie de quatre points 
source de reference (51) et on resout la formule de transformation (1), ce qui donne huit valeurs de 
coefficients de transformation a\ par serie de quatre points de reference (51) ; on considere le barycentre 45 
(54) de chacune des series de quatre points de reference (51) et on affecte a ce dernier les coefficients 
de transformation cti ainsi determines ; on considere des zones d'interpolation (T k ) de forme triangulaire 
delimitees par des segments de droite joignant lesdits barycentres (54) de series voisines, ces 
ban/centres constituant les sommets desdites zones triangulaires (T k ) ; on utilise la formule 
d 'interpolation suivante pour chaque coefficient de transformation a\ : 50 
(2) ai = pnCol + fr 2 Lig + p i3 

. on determine ainsi vingt-quatre coefficients d'interpolation <|3ij) par zone triangulaire d'interpolation (Tk) 
en resolvant la formule d'interpolation (2) pour chacun des coefficients de transformation (on) affecte a un 
sommet de zone triangulaire d'interpolation (Tk) ; 

ce procede consiste en outre, dans ladite phase de traitement, pour determiner les coordonnees 55 
spatiales (R, H) d'un point source dont on connaTt les coordonnees images (Col.Lig) de son image, a 
effectuer les operations suivantes : 

- rechercher, dans une operation de recherche, la zone d'interpolation (Tk) a laquelle le point appartient, 
ainsi que les coefficients d'interpolation (pjj) associes a cette zone, 

- calculer, dans une operation de caicul, les coordonnees spatiales (R, H) par application concomitante 60 
desdites formules de transformation (1) et d'interpolation (2). 

2. Procede selon la revendication 1 , caracterise en ce que lesdites operations de recherche et de caicul 
comportent les etapes suivantes, pour chaque point image traite : 

- une premiere etape d'initialisation (181), au cours de laquelle on attribue et on met en memoire des 
valeurs arbitraires (O, O) aux variables (R, H) que I'on cherche a determiner connaissant les coordonnees 65 
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images (Col,Lig) du point traite (Pi, n ). 

- une seconde etape d'initialisation (182), au cours de laquelle on initialise deux variables intermediaires 
(Ro, Ho) en mettant ces valeurs intermediaires egales a la derniere valeur connue des variables (R, H), on 
met a zero une variable k correspondant au numero de la zone d'interpolation (Tk), 

5 - une etape decrementation (183) au cours de laquelle on incremente la variable k, 

- une etape (184) de lecture des coefficients d'interpolation (Pij) associes a lazone d'interpolation T k , 

- un test d'appartenance (185) au cours duquel on recherche si le point ayant comme coordonnee 
lesdites variables intermediaires (R 0 , H 0 ) appartient a la zone d'interpolation portant le numero k (triangle 
Tk) ; dans la negative on incremente la valeur k et on recommence ce test ; 

10 - lorsque le test d'appartenance (185) est positif, au cours d'une etape de transformation (186), on vient 

lire dans un fichier image (Fla) les coordonnees images (Col,Lig) du point traite (Pi, n ), et on applique 
concomitamment les formuies de transformation (1) et d'interpolation (2) a partir des coefficients 
d'interpolation ((3ij) lus lorsdutest d'appartenance (185), pour calculer les coordonnees spatiales (R,H), 

- un test de conformite (187) au cours duquel on verifie I'egalite des coordonnees calculees a I'etape de 
15 transformation (186) avec les variables intermediaires (R 0l H 0 ) ; dans la negative, on met les variables 

intermediaires (R 0 , H 0 ) egales a celles calculees lors de I'etape de transformation et on revient a la 
seconde etape d'initialisation (182). 

3. Procede de reconstitution des coordonnees spatiales de chacun des points d'un ensemble de 
points source appartenant a une surface comportant une etape d'acquisition de donnees dans laquelle 

20 on eclaire, avec un faisceau plan de lumiere, la surface, on enregistre avec au moins une camera (13a, 

13b) au moins une image de la trace (p,p,9 n ) laissee par le faisceau plan sur la surface (20), la trace laissee 
par le faisceau plan de lumiere sur la surface etant dite "profit source", I'image de cette trace dans un plan 
image de la camera etant dite "profil image", on deplace relativement le faisceau plan et la surface et on 
enregistre avec ladite camera a nouveau au moins une image de la trace laissee par !e faisceau plan sur la 

25 surface, on recommence les operations de deplacement reiatif, d'eclairage et d'enregistrement en sorte 

que Ton enregistre ainsi une pluralite de profils Image constituant un echantillonnage de I'image de 
('ensemble de la surface, le procede etant caracterise en ce qu'il cornporte une etape de traitement au 
cours de laquelle, notamment, on considere les coordonnees lignes-colonne (Col, Lig) d'un ensemble de 
points image echantillonnant chaque profil image dans le plan image de ladite camera et les coordonnees 

30 spatiales des points source correspondant a chacun des points du plan image, cette etape de traitement 

comportant une phase d'homogeneisation au cours de laquelle, pour chacun des points traites (Pi, n ), on 
considere un nombre predetermine de points voisins, on calcule les coordonnees moyennes de ces 
points, et si le point traite est absent, ou si Tune au moins de ses coordonnees differe de la coordonnee 
moyenne correspondante d'un ecart superieur a une valeur predeterminee (E), on affecte au point traite 

35 lesdites coordonnees moyennes. 

4. Procede selon la revendication 3, caracterise en ce qu'au cours de la phase d'homogeneisation, pour 
chaque point source traite, dont on a determine les coordonnees spatiales (Ri, Hi, 0 n ), on effectue les 
operations suivantes : 

- on determine les points voisins du point (Rj, Hi, 8 n ) considere ; on considere les 2p + 1 profils source 
40 allant de 0 n -p a 0 n + P , 

- sur chaque profil source Gj ainsi considere, on selectionne les 2k + 1 points allant de (Ri-k, Hw<, 0j) a (Ri+k, 
Hi +k ,ei), 

- on calcule la moyenne R et la moyenne H des coordonnees (R, H) de chacun des points, 

- si Je point "traite" existe effectivement, on effectue une operation de filtrage au cours de laquelle pour 
45 chacun des points traites, on calcule une valeur absolue ( I Rj - R I ) ; si cette valeur est superieure a une 

valeur d'erreur donnee E, parametrable, on elimine ledit point traite (Ri, Hi, 0 n ), 

- si ie point (R,, Hi, B n ) "traite" n'existe pas ou s'il est elimine, on effectue un operation cMe de 
reconstitution au cours de laquelle on affecte au point (Ri, Hi, 0 n ), les valeurs moyennes calculees (R, H). 

5. Procede selon Tune quelconque des revendications 1 a 4, caracterise en ce que ledit axe directeur 
50 (Z) est un axe de rotation, tandis que ledit deplacement reiatif est une rotation autour de cet axe, une 

troisieme coordonnee spatiale (0) etant determinee par la coordonnee angulaire de la trace du faisceau 
de lumiere laisse dans un plan perpendiculaire (X, Y) audit axe directeur, cette troisieme coordonnee 
spatiale etant affectee a chaque profil enregistre par la camera (13a, 13b). 

6. Procede selon I'une quelconque des revendications 1, 2 et 5, caracterise en ce qu'il cornporte en 
55 outre une phase d'homogeneisation au cours de laquelle, pour chacun des points traites (Pi, n ), on 

considere un nombre predetermine de points voisins, on calcule les coordonnees moyennes de ces 
points, si le point traite est absent, ou si Tune au moins de ses coordonnees differe de la coordonnee 
moyenne correspondante d'un ecart superieur a une valeur predeterminee (E), on affecte au point traite 
lesdites coordonnees moyennes. 
SO 7. Procede selon la revendication 6, caracterise en ce que, au cours de la phase d'homogeneisation, 

pour chaque point source traite, dont on a determine les coordonnees spatiales (Ri, Hi ,9 n ), on effectue 
les operations suivantes : 

- on determine les points voisins du point (Ri, Hi, 0 n ) considere ; on considere les 2p + 1 profils source 
allant de0 n .p a 0n+ P) 

65 - sur chaque profil source Gj ainsi considere, on selectionne les 2k + 1 points allant de (Ri-k, Hi- k , 0j) a (Rj+k, 
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Hi+k.ej), 

- on calcule la moyenne R et la moyenne H des coordonnees (R, H) de chacun des points, 

- si le point "traite" existe effectivement, on effectue une operation de flltrage au cours de laquelle pour 
chacun des points traites, on calcule une valeur absolue I Ri - R I ; si cette valeur est superieure a une 
valeur d'erreur donnee E, parametrable, on elirnine ledit point traite (Ri, Hi, 8 n ), 5 

- si le point (R is Hi, 8 n ) "traite" n'existe pas ou s'il est elirnine, on effectue une operation dite de 
reconstitution au cours de laquelle on affecte au point (Ri, Hi, 8 n ), les valeurs moyennes calculees (R, H). 

8. Procede de realisation d'une image tridimensionnelle d'une surface tridimensionnelle mettant en 
oeuvre le procede selon Tune quelconque des revendications 1 a 7 applique a la determination des 
coordonnees spatiales d'un ensemble de points appartenant a ladite surface, caracterise en outre en ce 10 
qu'il met en oeuvre une machine-outi! adaptee a commander un outil seion cinq degres de liberte, a savoir 
trois degres de liberte permettant ie positionnement de cet outil dans I'espace, compte tenu desdites 
coordonnees spatiales et deux degres de liberte determinant un angle d'attaque de I'outil. 

9. Procede selon la revendication 8, caracterise en ce que ledit angle d'attaque est tel que I'outil est 
generalement dirige selon la normale a la surface image devant etre realisee et en ce que Ton procede a 15 
une detection des parties concaves de la surface pour corriger Tangle d'attaque de i'outil pour eviter 
autant que possible une interaction entre i'outil et la matiere usinee. 

10. Procede selon la revendication 7, caracterise en ce que pour la detection des parties concaves de la 
surface : 

- on considere la trace laissee par ladite surface (20) dans un plan perpendiculaire a un axe directeur (Z) 20 
et passant par le point en cours de traitement (Pi,n), cette trace etant dite secante (Si), 

- on considere egalement le profil source 0 n auquel appartient le point traite (Ri, Hi, 9n), 

- on detecte les parties concaves du profil en effectuant les operations suivantes : 

. on considere un vecteur normal ( DN^en.M, DN^enj+t) a chacun des deux points voisins (Rm, Hm, 8 n ) 

et (Ri+i, Hi+i, 9n), du point traite (Ri Hi, B n ) 25 

. on calcule le produit vectoriel desdits vecteurs normaux ( Dhf en.i-i, DN en,i+i) au profil aux points 
voisins du point traite ; on determine ainsi un coefficient de convexite QDj associe au point en traitement 

(Ri.Hi.Gn), 

. si le coefficient de convexite est negatif, le profil source est concave au point traite, et on calcule la 
somme (EiWi) desdits vecteurs normaux au profil, aux deux points voisins du point traite, 30 

( DN" 0n,i-1, 5N* 6nJ + 1 )' > 

. on memorise ce vecteur somme ( ed enj), 

. on considere alors I'ensemble des points traites appartenant a un meme profil (8 n ) ; on associe a chacun 
de ces points soit le vecteur normal au profil au point con sidere, lorsque le coefficient de convexite (QD) 
considere est positif, soit le vecteur somme, lorsque le coefficient de convexite (QD) est negatif, 35 
. lorsque I'ensemble des points appartenant a un meme profil est traite, on recommence un nombre 
predetermine de fois le traitement du profil en utilisant les vecteurs effectivement associes ( o^f e n ,i, 
""iU Bnj) au point traite (Ri, Hi, 8 n ). 

- on procede a un traitement analogue pour chacune des traces secantes (Si) et on associe a chacun des 
points traites (Ri, Hi, 8 n ) un vecteur correspondant soit a la normale a ladite trace secante Si au point 40 
traite ( SNen.O, soit a un vecteur somme, ( ES* enj) des vecteurs effectivement associes aux points 
voisins (Ri, H is 8 n -i) (R«, Hi, 8 n +i) au point traite (Ri, Hi, 8 n ). 

et pour chaque point, on dirige I'outil selon la direction des vecteurs effectivement associes audit point, 
dans le plan de son profil (8 n ) et dans celui de la trace secante (Si) a laquelle il appartient. 

11. Procede de realisation d'une image tridimensionnelle d'une surface tridimensionnelle dont on 45 
connait les coordonnees spatiales d'un nombre predetermine de points (Ri, Hi, 8 n ) echantillonnant cette 
surface, caracterise en ce qu'il met en oeuvre une machine-outil (18) adaptee a commander un outil selon 

cinq degres de liberte, a savoir trois degres de liberte permettant le positionnement de cet outil dans 
I'espace, compte tenu desdites coordonnees spatiales et deux degres de liberte determinant un angle 
d'attaque de I'outil de telle sorte que I'outil soit generalement dirige selon la normale a la surface devant 50 
§tre realisee et en ce que Ton procede a une detection des parties concaves de la surface, pour corriger 
I'angle d'attaque de I'outil pour eviter autant que possible une interaction entre I'outil et la matiere usinee, 
au cours de laquelle : 

- on considere la trace laissee par ladite surface (20) dans un plan perpendiculaire a un axe directeur (Z) 

et passant par le point en cours de traitement (Pi, n ), cette trace etant dite secante (Si), 55 

- on considere egalement un profil source 9 n auquel appartient le point traite (Ri, Hi, 8 n ), 

- on detecte les parties concaves du profil en effectuant les operations suivantes : 

. on considere un vecteur normal ( DNSnJ-i. DN*"eru+i) a chacun des deux points voisins (Ri-i,HM,8 n ) 
et (R i+ i, Hi+i,8n), du point traite (Ri,Hi,8n) _ 

. on calcule le produit vectoriel desdits vecteurs normaux ( nit en.M, DN en,i+i) au profil aux points 60 
voisins du point traite ; on determine ainsi un coefficient de convexite QDi associe au point en traitement 

(Ri, Hi8 n ), 

. si le coefficient de convexite est negatif, le profil source est concave au point traite, en on calcule la 
somme ( Iffi en.i) desdits vecteurs normaux au profil, aux deux points voisins du point traite, 
( DN* en.M, DN* 9n ' +1 )' 65 
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. on memorise ce vecteur somme ( ED^en.i), 

. on considere alors Pensemble des points traites appartenant a un meme profil (0 n ) ; on associe a chacun 
de ces points soit le vecteur normal au profil au point considere, lorsque le coefficient de convexite (QD) 
considere est positif, soit le vecteur somme, lorsque le coefficient de convexite (QD) est negatif, 
. lorsque I'ensemble des point appartenant a un meme profil est traite, on recommence un nombre 
predetermine de fois le traitement du profil en utilisant ies vecteurs effectivement associes ( l^f en.i, 

Er5" en.0 au point traite (R it Hi, 6 n ). 
- on procede a un traitement analogue pour chacune des traces secantes (Si) et on associe a chacun des 
points traites (R,-, Hi, 0 n ) un vecteur correspondant soit a la normale a ladite trace secante Si au point 
traite n ,i, soit a un vecteur somme, ( ESg n,«) des vecteurs effectivement associes (ES, SN) aux 

points voisins (R,, H it 0 n -i) (R., H h e n +i) au point traite (R (j H ( , G n ). 

et, pour chaque point, on dirige I'outil selon ia direction des vecteurs effectivement associes audit point, 
dans le plan de son profil (9 n ) et dans celui de la trace secante (Si) a laquelle il appartient. 
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